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ABSTRACT 

Thin-film  shape  memory  alloys  (SMAs)  have  become  excellent  candidates  for  microactuator  fabrication  in 
MEMS.  We  develop  a  material  model  based  on  a  combination  of  free  energy  principles  in  combination  with 
stochastic  homogenization  techniques.  In  the  first  step  of  the  development,  we  construct  free  energies  and 
develop  phase  fraction  and  thermal  evolution  laws  for  homogeneous,  single-crystal  SMAs.  Second,  we  extend  the 
single-crystal  model  to  accommodate  material  inhomogeneities  and  polycrystalline  compounds.  The  combined 
model  predicts  rate-dependent,  uniaxial  SMA  deformations  due  to  applied  stress  and  temperature.  Moreover, 
the  model  admits  a  low-order  formulation  that  is  suitable  for  subsequent  control  design.  We  illustrate  aspects 
of  the  model  through  comparison  with  thin-film  NiTi  superelastic  hysteresis  data. 

Keywords:  Shape  memory  alloy  model;  thin  film;  superelasticity;  activation  energy;  polycrystalline  com¬ 
pounds. 


1.  INTRODUCTION 

In  the  last  20  years,  shape  memory  alloys  (SMAs)  have  been  implemented  in  a  number  of  high  performance 
applications  requiring  high  work  densities,  large  recoverable  deformations,  and  high  stresses.  More  recently, 
shape  memory  alloys  have  become  excellent  candidates  for  microactuators  when  fabricated  as  thin  films.  Their 
work  output  density  (~50  MJ  m-3)  exceeds  that  of  other  microactuator  mechanisms  and  they  can  yield  large 
strokes  (~5%)  and  forces  (~30  mN).  Moreover,  thin-film  SMAs  heat  on  the  order  of  milliseconds  by  low  voltage 
(~5  V)  Joule  heating,  and,  unlike  bulk  SMAs  (wires,  beams,  etc.),  their  small  mass  and  large  surface-to-volume 
ratio  allow  fast  cooling,  potentially  permitting  switching  frequencies  on  the  order  of  100  Hz  [7]. 

Current  microscale  SMA  applications  include  microgrippers  [18],  micropumps  [29],  and  microcantilever 
switches  [16,  22].  Most  of  these  applications  rely  on  the  one-way  memory  effect  and  require  a  biasing  mechanism 
for  full  actuation.  However,  microdevices  using  functionally  graded  films  can  achieve  two-way,  out-of-plane 
displacement  with  a  smaller  footprint  than  conventional  micromechanisms  [11,  12].  Superelastic  NiTi  thin  films 
are  being  considered  for  high-strength  surface  coatings  in  MEMS  devices,  and  they  have  potential  for  microscale 
mechanical  energy  storage  devices,  and  vibration  dampeners  in  microelectronics  packaging  [13,  14].  The  reader 
is  referred  to  [19]  for  a  review  of  other  thin-film  SMA  applications. 

A  effort  to  develop  thin-film  SMA  material  models  is  warranted  to  ensure  the  existence  of  appropriate 
modeling  and  simulation  tools  needed  to  support  design,  performance  evaluation,  and  quantification  of  material 
capabilities.  Models  of  bulk  SMAs  capable  of  accommodating  nonlinear  material  behavior  pertinent  to  actuator 
design  are  available.  We  refer  the  reader  to  reviews  of  a  number  of  bulk  SMA  models  in  [4,  5,  6,  10,  25,  26]. 
In  conjunction  with  advancing  fabrication  and  film  characterization  technologies,  initial  thin-film  SMA  models 
have  been  developed  [2,  9,  17,  20].  Established  bulk  models  may  be  considered  for  thin-film  SMAs,  provided 
modifications  are  made  to  address  issues  of  scale  and  the  microfabrication  process  [24]. 

Focusing  on  thin-film  SMA  applications,  we  employ  as  a  starting  point,  the  Miiller-Achenbach-Seelecke 
theory  [1,  23,  27],  based  on  the  quantification  of  thermally  activated  processes  for  bulk  SMAs.  In  the  first 
step,  we  establish  local  free  energies  for  single-crystal,  homogeneous  SMAs  using  Boltzmann  statistics  to  derive 
rate  laws  for  phase  fraction  evolution.  In  addition,  we  formulate  a  balance  of  internal  energy  that  quantifies 
a  rate-dependent  release  of  latent  heats  and  heat  transfer  to  and  from  the  environment.  In  the  second  step 
of  the  development,  we  extend  the  single-crystal  model  to  accommodate  inhomogeneous  and  polycrystalline 
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materials  by  considering  material  parameters  to  be  manifestations  of  underlying  stochastic  densities.  The 
result  is  a  rate-dependent,  thermomechanical  model  that  predicts  relative  elongation  due  to  applied  stress  and 
temperature.  The  model  accommodates  SMA  behavior  pertinent  to  actuator  design  such  as  superelasticity  and 
the  shape  memory  effect.  Moreover,  the  model  admits  a  low-order  formulation  suitable  for  subsequent  control 
design,  and  most  of  the  model  parameters  are  identifiable  directly  from  standard  measurements.  In  Section 
2,  we  develop  the  homogeneous  SMA  model  and  we  extend  the  model  to  accommodate  nonhomogeneous  and 
polycrystalline  SMAs  in  Section  3.  Finally,  in  Section  4  we  validate  aspects  of  the  model  through  comparison 
with  superelastic  thin-film  SMA  data. 


2.  HOMOGENEOUS  MATERIAL  MODEL 

Motivated  by  the  theory  in  [27],  we  treat  a  lattice  layer  as  the  fundamental  element  of  our  model.  Considering 
deformations  along  one  axis  (across  a  single  habit  plane),  an  SMA  lattice  layer  admits  either  the  austenite  phase 
or  one  of  two  martensite  variants.  We  denote  the  volume  fraction  of  austenite  and  martensite  layers  in  an 
SMA  as  xa  (i),  where  a  represents  austenite  (A)  and  martensite  (M±).  The  phase  fractions  then  satisfy  the 
conservation  law 

Y^xa(t)  =  i  (i) 

a. 

over  all  time  t  >  0.  Throughout,  we  assume  the  martensite  variants  share  the  same  thermomechanical  proper¬ 
ties,  which  generally  differ  from  those  of  austenite.  Once  we  construct  a  thermoelastic  free  energy  relation  for 
SMAs,  we  will  model  the  evolution  of  the  phase  fractions  due  to  external  stresses  and  temperatures. 


2.1.  Free  Energy 


For  a  scalar  strain  e,  we  consider  the  potential  for  the  martensite  variants 

^(£iT)=^(£T£t)2  +  /1m(T),  (2) 

while  for  austenite  we  consider 

4A(e,T)=E^e2  +  (3A(T).  (3) 

The  strain-dependent  portions  of  the  potentials  represent  the  elastic  energies,  where  the  constants  Em  and  EA 
are  the  linear  elastic  moduli  for  martensite  and  austenite,  respectively.  The  quantity  St  corresponds  to  the 
stress-free  equilibrium  strain  of  martensite,  while  e  =  0  (no  deformation)  is  the  equilibrium  strain  for  austenite. 
The  parameter-dependent  family  of  functions 

Pa(T)  =  ca(T -TR) +ua-Tsa  (4) 


represent  the  chemical  (non-elastic)  free  energies  [9,  20,  27].  In  (4),  ua  denote  internal  energy  constants  and 
ca  are  phase-dependent  specific  heat  capacities  (constant  volume)  multiplied  by  the  material  density.  The  sa 
are  specific  entropies  of  the  form 


+  Vc 


(5) 


where  r]a  are  phase-dependent  entropy  constants.  In  (4)  and  (5),  Tr  is  the  temperature  of  the  reference  state 
from  which  energies  and  entropies  are  calculated. 

While  others  such  as  [9,  17,  20]  have  formulated  the  specific  Helmholtz  free  energies  as  a  mixture  of  (2)  and 
(3),  we  construct  a  single,  C1-continuous  Helmholtz  free  energy  by  joining  the  individual  potentials.  Specifically, 
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Figure  1.  Piecewise  quadratic  Helmholtz  free  energy  density  (6)  for  increasing  temperatures;  dashed  segments 
represent  the  concave,  unstable  states. 


we  employ  the  Helmholtz  relation 


i/j(e,T)={ 


^(£  +  £t)2 


~^P(£  +  £ 0  (T))2  +  Ipo  (T) 
^e2  +  A/3  (T) 
P(£-£o  (T))2  +  V>o(T) 

^(£-£Tf 


£  <  -JM  ( T ) 

-7m  (T)  <  e  <  -7a  (T) 
M  <  7a  (T) 

7 A  (T)  <  e  <  7m  (T) 
e  >  7 m  (T) 


(6) 


where  7m  (T)  and  7a  (T)  are  temperature-dependent  nodes  connecting  the  concave  parabolae,  which  represent 
unstable  states,  to  the  convex  potentials,  which  represent  stable  martensite  and  austenite  states.  To  facilitate 
subsequent  calculations,  we  have  shifted  the  martensite  minima  to  zero  so  that  the  austenite  minimum  has 
height 

A/3  (T)  =  pA  (T)  —  Pm  (T) .  (7) 

The  temperature-dependent  coefficients  Eg,  £g,  and  ip0  define  the  concave  parabolae  whose  maxima  are  ipg  at 
e  =  ±£q.  Enforcing  continuity  at  the  nodes  yields 


—73a 7 a  (T)  7M  (T)  -  Em  (7 m  ( T )  -  eT)  ( £t  -  7 a  (T))  =  2A/3  (T) ,  (8) 

which  implicitly  provides  a  thermodynamics-based  relation  for  the  nodes  via  (4)  and  (7). 

Note  that  the  local  three-well  energy  (6)  is  only  valid  at  temperatures  where  austenite  is  stable.  We  define 
a  local  critical  temperature  Tc  below  which  the  austenite  phase  is  unstable  upon  cooling.  Then,  (6)  holds  only 
for  T  >  Tc.  At  the  critical  temperature,  the  austenite  nodes  converge  to  7a  (Tc)  =  0,  thereby  eliminating 
the  austenite  potential  well  for  all  lower  temperatures.  The  Helmholtz  free  energy  reduces  to  a  double-well 
potential  for  T  <  Tc.  In  this  paper,  we  focus  on  model  development  for  the  three-well  potential.  Figure  1 
illustrates  (6)  for  the  full  range  of  temperatures  under  the  assumption  =  sm  —  sa  <  0. 

The  Gibbs  free  energy  density  associated  with  the  Helmholtz  free  energy  is 

G  (e,  cr,  T)  =  i))  (e,  T)  -  ere,  (9) 
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Figure  2.  The  specific  Gibbs  free  energy  (9)  at  a  fixed  temperature  and  varying  stress.  The  equilibrium 
conditions  (10)  yield  a  hysteretic  stress-strain  relation  with  jumps  at  the  nodes. 


where  a  is  the  stress  conjugate  to  e.  Given  a  stress  a  at  temperature  T,  stable  equilibrium  strains  must  satisfy 


dG  n  /G  n 
—  =  0  and  — -y-  >  0. 
oe  oe z 


(10) 


Given  the  quadratic  form  of  the  local  elastic  energy,  the  equilibrium  stress-strain  relationship  yields  a  linear 
Hooke’s  law  with  hysteresis,  as  illustrated  in  Figure  2.  The  strain  values  at  which  jumps  occur  coincide  with 
the  nodes  ±7m  (T)  and  ±7,4  (T),  since  they  are  also  inflection  points  of  the  Helmholtz  energy.  We  define 


a  a  ( T )  =  Ea'Ja  ( T ) 


(11) 


and 

ctm  (T)  =  Em  (7 m  ( T )  —  £T) , 

corresponding  to  the  stress  values  at  the  jumps.  As  in  [27],  we  assume  the  local  relation 

(T)  -  aM  (T)  =  S , 


(12) 

(13) 


where  6  is  a  material  dependent  constant  estimated  from  hysteresis  data.  By  the  construction  of  the  Gibbs 
potential,  a  a  (T)  >  &m  (T)  for  all  T,  so  <5  >  0.  With  (11)  and  (12),  we  use  (13)  to  formulate  (8)  in  terms  of  a 
single  temperature-dependent  stress  (or  strain).  For  example, 


A/3  (T) 


6£T  2Em5£T  -  (1  -  Er)  [5  +  (T)] 

c\  O  A  (  )  r>  7 


(14) 


where  Er  =  Em/Ea  is  less  than  one  for  shape  memory  compounds.  Note  that,  since  a  a  (Tc)  =  0,  (14)  provides 
an  implicit  expression  for  the  critical  temperature  solely  in  terms  of  material  constants. 

The  first  order  Taylor  expansion  of  (7)  about  T  =  Tr  is 


A/3  (T)  =  A u-  ArjT 


(15) 


with  truncation  error  O  yT  —  Tr)2J  ,  where  A  u  =  ua  —  um  and  A77  =  r)A  —  t]m- 

relation  when  T  =  Tr  and  for  materials  in  which  Ac  =  ca  —  cm  =  0.  Using  (14) 
expression  for  Tc: 

T  ~  -  fer 

c  Arj  2Arj 


Note  that  (15)  is  an  exact 
and  (15),  we  get  an  explicit 


(16) 
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In  an  SMA,  individual  austenite  layers  become  unstable  upon  cooling  at  temperatures  between  the  measured 
martensite  start  Ms  and  martensite  finish  Mf  transformation  temperatures.  Therefore,  we  estimate  Tc  for  a 


layer  using  the  relation 


Ms  +  Mf 
2 


(17) 


With  measurements  of  transformation  temperatures  and  a  transformation  stress  (e.g.,  obtained  via  a  tensile 
test),  we  use  (14)  and  (17)  to  calculate  the  chemical  free  energy  constants  Aw  and  A for  an  alloy.  For  example, 


[2EM£t  -  (1  -  Er)  8\  ah  +  {  1  ~  Er)  a2h 
2Em  (Th  -  Tc) 

Aw  =  AijTc  + 


where  Th  >  Tc  is  the  temperature  at  which  loading  transformation  stress  ah  =  a  a  ( Th )  is  measured.  A  similar 
approach  for  identifying  the  energy  constants  is  employed  in  [20].  Ultimately,  we  are  able  to  define  the  Gibbs 
potential  completely  from  material  data. 


2.2.  Phase  Evolution 

Based  on  the  local  free  energy,  we  follow  the  approach  in  [27]  and  model  the  phase  evolution  with  the  rate  laws 

x_  (t)  =  Pa-Xa  (i)  -  P-X-  (i) 

xA  (t)  =  P-X-  ( t )  -  Pa-xa  ( t )  +  P+x+  (t)  -  Pa+xa  (t)  (20) 

x+  ( t )  =  Pa+xa  (t)  -  P+x+  (t) , 

where  *  denotes  differentiation  in  time.  The  functions  P±  denote  the  likelihoods  that  lattice  elements  will 
undergo  a  transformation  (to  either  austenite  or  a  martensite  variant),  while  Pa±  denote  the  likelihoods  that 
austenite  will  transform  to  M  .  Using  the  conservation  relation  (1),  the  rate  law  reduces  to  the  two  coupled 
ODEs 

x-  ( t )  =  -  (P-  +  Pa-)  x-  (t)  -  Pa-x+  ( t )  +  Pa - 

x+  (t)  =  -  (P+  +  Pa+ )  a;+  (t)  -  Pa+X-  (t)  +  PA+ ■  (21) 

Next,  we  formulate  the  likelihood  functions  used  in  (21). 


2.3.  Transformation  Likelihoods 

We  quantify  the  probability  that  a  layer  will  attain  specific  energy  G  using  the  Boltzmann  relation 

H{G)  =  Ce~GV/kBT,  (22) 

where  C  is  a  normalization  factor  chosen  to  yield  a  probability  of  one  when  integrating  n  (G)  over  all  energy 
states.  In  (22),  V  is  the  volume  of  the  layer,  and  ks  is  the  Boltzmann  constant.  Using  (22),  we  compute  the 
likelihood  of  a  martensite  layer  undergoing  a  transformation. 

P±(a,T) 

given  a  stress  a  and  temperature  T,  where  m  is  the  mass  of  the  layer  and  G  (±7 M,  <7,  T)  is  the  Gibbs  energy  (9) 
at  nodes  ±7 m  (T).  We  have  set  the  integration  limits  for  the  normalization  factor  to  cover  all  possible  stable 
(or  metastable)  martensite  equilibrium  states,  neglecting  the  unstable  states  defined  by  the  concave  parabolae 


0-G(±7M,<y,T)V/kBT 


2tt mV2/3  r±0°  e-GG,<r,T)V/kBTd~ 
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in  (6).  The  factor  in  front  of  the  normalized  Boltzmann  probability  represents  the  frequency  of  attempts  a 
layer  makes  to  overcome  an  energy  barrier.  Naturally,  the  frequency  increases  with  increasing  activation  energy, 
quantified  by  the  product  kET .  Likewise,  for  the  likelihood  of  an  austenite  layer  undergoing  a  transformation, 
we  have 


Pa±  (a ,T) 


kBT 


0-G(±7A,<y,T)V/kBT 


2irmV2/3  f™A  e-G(7^T)v/feBTd7' 


(24) 


Using  (11),  (12),  and  (14),  we  simplify  the  normalization  integrals  and  express  the  likelihoods  in  terms  of 
transformation  stresses. 


P±{cr,T) 


\J  Er 


g— (crM=FCT)2/2<^M(r)2 


t  erfc  [(ctm  T  cr)  /a/2 toM  (T)] 


(25) 


and 


Pa±  (c,  T) 


1  g-(<7ATU2  /2ua(T)2 

t  erf  [(a a  +  cr)  / a/2 loa  (T)]  +  erf  [{aA  -  a)  / y/2u>A  (T)] 


where 


is  the  standard  error  function  and  erfc  (x) 


erf  ( x )  =  —j= 


dr 


=  1  —  erf  (x)  is  its  complement. 


The  resulting  factor 


(26) 

(27) 


r  (28) 
has  units  of  seconds  and  represents  a  transformation  relaxation  time.  The  functions 

ua(T)  =  ^J^  (29) 

quantify  thermal  activation  energy  densities  for  martensite  and  austenite.  Both  (25)  and  (26)  have  qualities 
similar  to  those  of  normal  distributions,  where  the  transformation  stresses  aa  ( T )  act  as  means  and  uia  ( T )  as 
standard  deviations.  Increasing  activation  energy  densities  (large  variances)  allow  transformations  to  occur 
at  wider  ranges  of  stress  values  approaching  aa  (T).  Diminishing  activation  energy  densities  (small  variances) 
only  allow  transformations  to  occur  at  stresses  very  near  aa  (T). 


2.4.  Thermal  Evolution 


We  describe  uniform  temperature  changes  in  the  material  via  a  simplified  balance  of  internal  energy. 

c  (t)  T  (t)  =  -hcQ  [T  -  Te  (t)]  -  ^2  haxa  +  (30) 

ol 


where  c (t)  =  ^2acaxa  ( t )  represents  the  average  phase-dependent  specific  heat.  The  right-hand  side  of  (30) 
accounts  for  heat  convection  to  the  environment,  where  hc  is  a  convection  coefficient,  U  is  the  SMA  surface 
area-to- volume  ratio,  and  TE  (t)  is  the  external  temperature.  In  general,  hc  is  temperature  and  geometry 
dependent.  For  this  paper,  we  assume  hc  is  constant  and  we  refer  the  reader  to  [3,  7,  24]  for  discussions  of 
variable  hc  pertinent  to  bulk  and  thin-film  SMAs. 

Equation  (30)  also  accounts  for  rate-dependent  heat  generation  and  loss  due  to  phase  transformations,  where 
the  specific  enthalpies  ha  have  the  form 

ha=  ga  +  Tsa.  (31) 

In  (31),  ga  is  a  local  minimum  of  (9)  and  sa  is  the  specific  entropy  from  (5).  Given  stress  cr  at  temperature  T, 


9A  =  dr  +  A/3(T) 

(32) 

-a2  .  . 

=  T G£T * 

(33) 
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Since  ]Pa  xa  ( t )  =  0,  we  have 


(34) 


haxa  =  ( h _  —  Ha)  X-  +  (h+  —  h a)  x+, 

oc 

and 

h±  -  hA  =  ~4 T—  ( vr - T  o-ffT  -  Ac (T  -  Tr)  -  Au.  (35) 

A  V  “AJ  / 

The  differences  in  martensite  and  austenite  specific  enthalpies  in  (35)  are  referred  to  as  latent  heats  of  trans¬ 
formation. 

Lastly,  in  (30),  J  (t)  incorporates  heat  generation  via  resistive  Joule  heating,  where 

J  (t)  =  pe  ( t )  (36) 

for  the  average,  phase-dependent  electrical  resistivity  pe  (t)  =  Paxa  {t)>  applied  electric  current  /,  and 
the  SMA  cross  sectional  area 


2.5.  Average  Strain 


Given  the  phase  fractions  and  the  thermal  evolution  modeled  by  (21)  and  (30),  we  quantify  the  physical  response 
to  stress  and  temperature.  Boltzmann  statistics  govern  the  response  of  individual  martensite  and  austenite 
layers.  For  martensite,  the  expectation  strains  are 


r-lM 


{£-)  = 


7M(G(7,ct,  T^d'y 


(37) 


and 


7/r  (G(7,cr,  T))dj, 


(38) 


where  the  Boltzmann  function  p  (G)  from  (22)  acts  as  a  probability  density  function.  Likewise  for  austenite, 
we  have 

/'I  A 

IP  (G  (7,  cr,  T))  d'y.  (39) 

-7  A 

As  with  the  likelihood  functions,  we  have  set  the  integration  limits  to  cover  all  possible  stable  martensite 
equilibrium  states,  neglecting  the  unstable  states  defined  by  the  concave  regions  in  (6).  Evaluating  the  integrals 
and  using  (11),  (12),  and  (14)  yield 


(e±) 


±tujm  ( T ) 

EmVE  r7 r 


P±(v,T)  + 


(40) 


and 

(sA)  =  [Pa-  (a,  T)  -  PA+  (a,  T)]  +  ^  (41) 

in  terms  of  the  likelihood  functions  (25)  and  (26).  From  (40)  and  (41),  it  is  apparent  that  the  expectation 
strains  are  the  equilibrium  solutions  (10)  to  the  local  Gibbs  free  energy  with  perturbations  whose  magnitudes 
depend  on  the  activation  energy.  Indeed,  in  the  limit  of  vanishing  activation  energy  densities  u>a  (T),  (40)  and 
(41)  converge  to  the  equilibrium  solutions. 

With  the  expectation  strain  of  a  layer  in  each  phase,  we  calculate  the  average  local  strain  as  the  weighted 
sum 

£mech  —  y  )  Xa  {^a)  •  (4^) 

a. 
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Apart  from  the  strain  £mech  attributed  to  mechanical  deformation  and  phase  transformation  processes,  there 
are  strains  due  to  thermal  expansion  as  the  temperature  varies  over  time.  We  model  the  thermal  strain  as 

£therm  =  A  (t)  (T  —  Tq)  (43) 

for  the  average,  phase-dependent  thermal  expansion  coefficient  A  (t)  =  '}2a  Xaxa  (t),  and  initial  temperature 
To.  Therefore,  the  total  average  strain  of  a  layer  in  response  to  a  combination  of  applied  stress,  external 
temperature,  and  current  is 

£  —  £-mech  4“  £  therm-  (44) 

Typically,  £therm  is  negligible  compared  to  £mech  for  moderate  temperature  changes  as  would  be  encountered 
in  bulk  superelastic  applications.  However  for  thin  films  in  MEMS  applications,  thermal  strains  can  have 
significant  effects  on  microdevice  behavior  (e.g.,  see  [16,  22]). 

We  have  modeled  the  phase  fractions  and  expected  strains  with  free  energy  expressions  that  are  valid  for  local 
lattice  behavior.  Nevertheless,  the  model  does  quantify  macroscopic  strains  for  homogeneous,  single-crystal 
SMAs  in  which  the  bulk  lattice  exhibits  uniform  local  behavior.  In  the  next  section,  we  extend  our  uniform 
lattice  model  to  inhomogeneous  compounds,  where  the  local  thermomechanical  behavior  can  vary  throughout 
the  material. 


3.  INHOMOGENEOUS  MATERIAL  MODEL 


Most  thin-film  SMAs  are  polycrystalline  due  to  the  sputter  deposition  process.  Moreover,  precipitates  and 
process-induced  impurities  may  pervade  the  films.  In  general,  poly  crystalline  and  inhomogeneous  materials  have 
a  nonuniform  lattice  structure  that  can  yield  variations  in  the  local  material  response  as  well  as  a  nonuniform 
stress  field.  Measurements  of  macroscopic  phenomena,  such  as  stress-strain  hysteresis  and  transformation 
temperatures,  reflect  an  average  of  these  variations.  Accordingly,  to  account  for  inhomogeneities,  we  treat 
measured  quantities,  such  as  £t  and  5,  as  manifestations  of  distributions,  rather  than  fixed  values.  In  this 
paper,  we  construct  a  statistical  distribution  about  5  by  considering 

poo 

£(cr,T)=/  £(<j,T;5)  fi(6)d5:  (45) 

Jo 

where  £  is  the  local  strain  in  (44),  and  fi  is  a  probability  density  function  defined  on  [0,  oo).  Note  that  variations 
in  the  local  transition  temperature  Tc  are  also  accommodated  by  varying  6  via  (14)  and  (16). 

To  include  the  effects  of  a  nonuniform  stress  field  associated  with  an  inhomogeneous  lattice,  we  represent 
the  effective  stress  on  a  sample  as  a  perturbation  about  the  macroscopic  applied  stress.  In  this  case,  we  take 


e(cr,  T) 


£(ae,T)  f2(a  -  ae)  dae 


(46) 


where  ae  is  the  effective  stress  perturbation  and  f2  is  a  probability  density  function  defined  on  (—00,00).  In 
particular,  by  distributing  about  the  applied  stress,  we  account  for  polycrystal  grains  undergoing  stress-induced 
transformations  at  different  stress  levels.  We  note  that  a  trigonometric  distributional  averaging  is  proposed 
in  [23]  to  accommodate  stress  variations  on  rotated  polycrystal  grains.  The  combination  of  (45)  and  the 
convolution  of  (46)  yields  a  relation  for  the  total  macroscopic  response  to  a  time- varying  stress  and  an  evolving 
internal  temperature 

pOO  poo 

£(<7,T)=  /  £(a  -  ae,T;5)  fi(6)  f2((je)  daed6. 

J  0  J  —00 

In  the  next  section,  we  validate  the  full  model  (47)  with  thin-film  SMA  hysteresis  data. 


(47) 


4.  MODEL  VALIDATION 


A  NiTi  film  of  8/iin  thickness  was  deposited  onto  a  Si  substrate  using  a  dedicated  DC  magnetron  sputtering  sys¬ 
tem  built  at  the  UCLA  Active  Materials  Lab  [28].  As  described  in  [12],  the  system  has  a  temperature-controlled 
target  heater  used  to  make  functionally  graded  films.  The  target  was  rated  47.0at%Ni,  but  microprobe  (JEOL 
JXA-8200)  tests  prior  to  deposition  indicated  the  target  was  48.5at%Ni  and  the  finished  polycrystalline  film 
measured  48.8at.%Ni.  After  the  film  was  peeled  from  the  substrate,  it  was  post-annealed  at  500  °C  for  120 
minutes. 

After  annealing,  the  transformation  temperatures  were  measured  with  a  Shimadzu  differential  scanning 
calorimeter  yielding  As  =  56  °C,  Af  =  80  °C,  Ms  =  61  °C,  Mf  =  50  °C.  Tensile  specimens  having  the 
dimensions  0.1778  cm  width  and  1.5  cm  gauge  length  were  cut  from  the  film.  Using  an  MTS  microtensile 
apparatus,  load-controlled  (0.94  MPas-1)  tensile  tests  were  performed  in  a  thermal  chamber.  The  martensite 
modulus  Em  =  14.0  GPa  was  measured  from  an  initial  tensile  test  at  25  °C.  The  results  for  the  full  cycle 
test  conducted  at  80  °C  are  plotted  in  Figure  3.  Near-superelasticity  is  observed  with  a  0.11%  residual  strain. 
The  austenite  modulus  Ea  =  22.4  GPa  was  extrapolated  from  the  loading  curve,  and  the  measured  loading 
transformation  ah  =  195  MPa. 

From  the  superelastic  curve  at  80  °C,  we  estimated  et  =  0.0233  by  extrapolating  a  tangent  line  from  the 
unloading  curve  and  estimated  the  mean  5  =  125  MPa  from  the  hysteresis  thickness.  We  did  not  measure 
the  specific  heats  and  published  values  for  thin-film  SMAs  were  not  available  as  estimates.  As  is  typical 
for  bulk  NiTi,  we  assume  Ac  ss  0.  Therefore,  (30)  and  (35)  simplify  accordingly,  and  ah,  Ms,  and  Mf  are 
sufficient  to  estimate  (7).  We  calculated  hc  =  140  Wm-2K-1  from  [24],  and  the  surface  area-to-volume  ratio 
=  2.511  x  105  m-1.  We  computed  a  log-normal  distribution  about  5  using  a  variance  of  (27  MPa)  ,  and 
we  used  an  effective  stress  variance  of  (5  MPa)  with  a  normal  distribution.  For  the  average  layer  volume, 
we  took  V  =  5.0  x  10-23  m3  and  chose  a  relaxation  time  r  =  1.7  ms.  For  the  remaining  parameters,  we 
used  published  NiTi  values:  p  =  6450  kgm-3,  cVA  =  cVM  =  837.3  Jkg-1K-1,  A^  =  11.0  x  10-6  K-1,  and 
Am  =  6.6  x  10-6  K-1.  For  alternative  approaches  to  identifying  and  measuring  these  parameters,  refer  to 


Figure  3:  Model  fit  to  thin-film  NiTi  superelastic  hysteresis  data. 
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[3,  8,  21].  The  results  of  the  full  model  with  the  identified  parameters  is  compared  with  the  data  in  Figure  3. 
It  is  observed  that  the  model  accurately  predicts  the  asymmetric  loading-unloading  curves. 


5.  CONCLUDING  REMARKS 

The  theory  presented  here  provides  a  technique  for  modeling  uniaxial  hysteresis  in  SMAs.  In  the  first  step 
of  the  development,  we  employ  the  theory  of  Miiller-Achenbach-Seelecke  to  develop  a  local  rate-dependent, 
thermomechanical  model  applicable  to  homogeneous  single  crystals.  Then  we  accommodated  inhomogeneous 
and  polycrystalline  materials  by  averaging  the  variations  of  local  material  parameters.  The  full  macroscopic 
model  accurately  characterizes  superelastic  hysteresis  of  thin-film  SMAs  as  demonstrated  in  Figure  3.  Moreover, 
the  model  admits  a  low-order  formulation  that  is  suitable  for  subsequent  control  design.  Future  work  includes 
the  development  of  control  algorithms  employing  robust  and  multi-objective  control  frameworks.  This  will 
build  on  previous  control  analysis  for  piezoceramic  and  magnetostrictive  materials  employed  in  nonlinear  and 
hysteretic  regimes.  In  addition,  after  further  validation,  we  will  consider  integrating  the  material  model  into 
a  finite  element  analysis  infrastructure  to  support  future  design,  performance  evaluation,  and  qualification  of 
thin-film  SMA  microcomponents. 
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